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PLASTIC  STRESS -STRAIN  RELATIONS  BASED  ON  INFINITELY 
MANY  PLANE  LOADING  SURFACES* 
by 
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ABSTRACT 


This  paper  is  concerned  with  the  development  of  a 
theory  of  plastic  stress-strain  relations  for  work  hardening 
materials  based  on  infinitely  many  plane  loading  surfaces.  The 
stress-strain  relations  belonging  to  this  class  are  closely 
related  to  those  the  the  linear  incremental  type  but  have  the 
property  of  being  integrable  in  a restricted  sense.  They  are 
also  non-linear  in  that  a corner  appears  in  the  yield  surface 
at  the  point  of  loading.  A stress-strain  relation  of  this  type 
for  isotropic  materials  is  presented.  The  problem  of  including 
a description  of  Bauschinger  and  allied  effects  within  the 
theory  is  considered. 


The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  N7onr-35801  with  Brown  University. 
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Introduction 

Some  recent.  d ^velopments  [1],  [2]  have  made  it  feasible 
to  investigate  the  possibilities  of  a certain  class  of  non- 
linear incremental  stress-strain  relations  for  work  hardening 
materials.  Linear  incremental  stress-strain  relations  have  been 
extensively  investigated  and  results  of  a quite  general  nature 
have  been  achieved  [3  ],  [ *+  ] . Studies  of  the  behavior  of  metals 
from  a phenomenological  point  of  view  ha\?e  revealed  a few  condi- 
tions which  greatly  restrict  the  possible  forms  of  stress-strain 
relations,  linear  or  not.  In  particular  a work  hardening  condi- 
tion postulated  by  Drucker  [2]  has  for  reaching  consequences. 

Under  this  condition  he  has  shown  that  the  yield  surface  must  be 
convex,  that  at  any  smooth  point  of  the  yield  surface  the  strain 
increment  vector  must  be  normal  to  the  surface  for  loading  from 
that  point,  and  that  there  are  certain  restrictions  on  the  direc- 
tions of  the  strain  increment  vector  for  loading  from  a corner. 
From  these  last  restrictions  it  follows  that  the  stress-strain 
relations  must  be  non-linear  for  loading  from  a corner. 

Before  these  latter  results  were  obtained  a stress- 
strain  relation  based  on  physical  consideration  of  a crystal 
aggregrate  had  been  proposed  by  Batdorf  and  Budiansky  [ 5 ] . It 
was  later  shown  to  satisfy  Drucker' s worv  hardening  condition 
and  to  possess  a uniqueness  and  a variational  theorem  [ 6 ] . This 


* 


Numbers  in  square  brackets  refer  to  Bibliography  at  the  end 
of  the  paper. 
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theory,  called  the  slip  theory  of  plasticity,  is  non-linear  in 
that  a corner  is  always  formed  in  tho  yield  surface  at  the  load- 
ing point.  Since  the  stress-strain  relations  of  slip  theory 
were  formulated  without  mention  of  a yield  surface  and  since 
practically  all  previous  theories  had  been  linear  (and  hence 
properly  applicable  only  in  case  of  a smooth  yield  surface)  it 
v/as  thought  that  slip  theory  stood  in  a class  by  itself. 

The  wor1,  done  to  obtain  theorems  regarding  the  nor- 
mality of  the  strain  increment  vector  to  the  yield  surface  was 
partly  motivated  by  a desire  to  use  the  corresponding  yield 
function  (or  loading  function)  as  a plastic  potential  in  the 
stress-strain  relations.  Recentlv  :oiter  [l  ] showed  that  the 
concept  of  plastic  potential  may  be  retained  even  if  the  yield 
surface  is  singular  (with  corners,  ridges,  etc.)  by  the  simple 
device  of  introducing  more  than  one  loading  function.  The  yield 
surface  is  the  boundary  of  those  points  in  stress  space  which 
represent  clastic  behavior  with  respect  to  all  loading  surfaces 
and  of  course  the  yield  surface  is  singular  whore  tv/o  or  more 
loading  surfaces  intersect.  In  such  a theory  of  plasticity  the 
total  plastic  strain  increment  ig  the  sum  of  the  contributions 
from  each  loading  function.  The  stress-strain  relation  is 
necessarily  non-linear  at  corners  in  the  yield  surface  but  may 
possibly  be  non-linear  only  at  the  corners.  Certain  uniqueness 
and  variational  theorems  were  extended  by  Koiter  to  cover  this 
class  of  non-linear  theories.*  He  also  showed  that  slip  theory 

Which  also  satisfy  the  work  hardening  condition  and  the  so- 

called  conditions  of  continuity  and  consistency  [3  ] . 
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belongs  to  this  class.  The  loading  surfaces  of  slip  theory  are 
an  infinite  set  of  planes. 

A qualitative  discussion  of  theories  of  plasticity 
based  on  many  plane  loading  surfaces  was  given  in  [7  ].  It  was 
shown  that  in  certain  cases  the  resulting  stress-strain  rela- 
tions are  integrable  in  a restricted  sense  and  thus  partially 
resemble  the  stress-strain  relations  of  deformation  theories  of 
plasticity.  The  present  paper  is  concerned  with  further  develop- 
ment of  the  theory.  A stress-strain  relation  for  initially 
isotropic  materials  is  constructed  and  certain  interesting 
special  cases  are  discussed.  It  is  shown  how  an  account  of 
Bauschinger  and  allied  effects  may  be  included  within  the  theory 
by  postulating  a relation  between  the  displacements  of  the  in- 
finitely many  loading  planes  which  envelope  the  yield  surface. 

A function  depending  on  loading  history  alone  is  introduced 
which  determines  both  the  yield  surface  and  the  plastic  strain. 

Plane  Loading  Surfaces  ir.  General . 

We  begin  by  noting  that  at  least  some  planes  in  stress 
space  have  a direct  physical  interpretation.  The  shear  stress 
in  the  direction  X^  on  an  element  of  area  normal  to  is  given 
by 


x = dij  X}  pj  (1) 

where  X^  and  are  unit  vectors.  The  geometrical  representation 
in  stress  space  of  the  equation  x - constant  is  of  course  a plane. 
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Suppose  that  the  critical  shear  stress  on  a certain  slip  system 
in  a single  crystal  is  k.  In  this  case  % - k might  be  regarded 
as  the  initial  yield  condition  for  the  slip  system,  the  corres- 
ponding plane  in  stress  space  as  the  initial  yield  surface,  and 
the  function  t of  the  stresses  (1)  as  the  loading  function.  One 
may  easily  verify  that  the  shear  strain  due  to  slipping  is  repre- 
sented by  a vector  normal  to  this  plane  (in  nine  dimensional 
stress  space).  If  a plane  loading  surface  is  introduced  for  all 
possible  orientations  of  slip  systems  in  a random  aggregate  of 
crystals  and  certain  othe^  assumptions  are  made  than  the  result 
is  the  so-called  slip  theory  of  plasticity  [ 5 ] • In  this  paper 
the  theory  of  stress-strain  relations  for  work  hardening  materials 
based  on  an  infinite  number  of  plane  loading  surfaces  is  con- 
sidered as  a.  subject  in  itself.  No  attempt  will  be  made  to 
develop  a rational  theory  based  on  the  physics  of  crystal  aggre- 
gates but  some  suggestions  from  that  source  will  be  of  use. 

Suppose  that  an  initial  yield  surface  for  a strain 
hardening  material  has  been  given.  Such  a surface  must  be  convex 
hence  it  can  be  regarded  as  an  envelope  of  planes  which  nowhere 
intersect  the  surface.  Tror  the  present  assume  that  each  plane 
behaves  as  though  the  others  wore  not  there.  As  loading  into  the 
plastic  range  proceeds  each  nlnne  is  moved  outward  parallel  to 
itself  by  the  loading  point  more  or  less  or  not  at  all  depending 
on  the  loading  path.  Two  examples  of  the  resulting  yield  bound- 
ary for  biaxial  loading  are  shown  in  Figs.  1(a)  and  1(b).  Each 
plane  remains  as  far  out  as  it  has  ever  been  pushed  by  the 
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loading  point  during  the  course  of  loading  and  ordinarily  there 
is  a corner  in  the  yield  surface  at  the  final  load  point.*  At 
each  stage  of  loading  the  yield  surface  is  convex,  as  it  must  be 
for  a work  hardening  material  according  to  a theorem  of  Drucker 
[ 2 ]. 

The  plastic  strain  is  assumed  to  be  the  sum  of  con- 
tributions from  each  loading  plane.  The  contribution  from  a 
single  plane  is  assumed  to  be  given  by  a linear  stress-strain 
relation  of  the  Prager -Drucker  type  in  which  the  loading  function 
f is  linear  and  homogeneous  in  the  stresses  and  where  f = c > 0 
is  the  equation  of  the  loading  plane.  The  contribution  of  a 
group  of  planes  whose  normals  lie  within  a snail  solid  angle  is 
given  by: 

b = G(f,w)  n.j  f f > 0 

(2) 

= 0 f < 0 

rVF 

where  n^  = ~ is  a constant  normal  vector  to  one  of  the  load- 
13  a°ij 

ing  planes,  6u  is  an  element  of  solid  angle,  and  G > 0 if  f 
exceeds  a certain  value  which  is  the  yield  point  for  this  plane 
but  vanishes  otherwise.  The  possible  dependence  of  G upon  the 
orientation  of  the  loading  plane  is  indicated  by  writing  it  as 
a function  of  u . Equation  (2)  nay  be  integrated  with  respect 
to  time  to  give: 


* 


See  [ 8 ] for  recent  experimental  evidence  of  corners 
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6 EjJ  = H(f,w)  n^  6w  (3) 

(where  -SH  = G)  as  may  easily  be  verified  by  differentiating  (3) • 
0f 

The  total  plastic  strain  is  obtained  by  integrating  over  all 
orientations  of  loading  planes: 


e 


P 

ij 


(4) 


The  value  of  f apnearing  in  (3)  is  the  largest  value  attained 
during  loading  and  is  a measure  of  the  distance  the  loading 
plane  has  been  pushed  out.  As  a result  of  the  integrability  of 
(2)  the  strain  contributed  by  each  loading  plane  depends  on  the 
loading  path  only  in  so  far  as  it  depends  on  the  distance  the 
loading  plane  has  been  pushed  out  during  the  course  of  loading. 
As  a consequence  of  this  integrability  the  following  statement 
may  be  made  concerning  the  stress-strain  relation  (4):  l,The 

total  plastic  strain  is  the  same  for  any  two  loading  paths  which 
result  in  the  same  yield  surface.”  In  this  sense  the  stress- 
strain  relation  is  path  independent.*  This  of  course  is  not  the 
same  as  saying  that  the  total  plastic  strain  is  the  same  for  any 
two  loading  paths  which  reach  the  same  final  load  point. 

Any  loading  path  resulting  in  a yield  surface  which 
could  also  be  produced  by  a radial  loading  path  will  be  called 
a "nearly  radial  loading  path."  In  Fig.  2 the  nearly  radial 
loading  path  OPQ  results  in  the  same  yield  surface  as  the  radial 
path  OF'Q.  Note  that  a path  which  is  not  nearly  radial  (OP"P' ) 


* 


In  the  case  of  slip  theory  this  theorem  was  known  to  its  authors 
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may  possibly  be  made  so  by  further  loading  (OP"?^).  In  Fig.  3 
the  loading  point  has  reached  P by  a nearly  radial  loading  path. 

If  the  loading  path  proceeds  into  the  shaded  region  then  all 
planes  previously  loaded  continue  to  load.  This  region  -will  be 
called  the  region  of  total  loading.  Any  loading  path  which 
always  proceeds  into  the  region  of  total  loading  is  nearly  radial. 

So  far  it  has  been  assumed  that  each  plane  acts  in- 
dependently of  the  others  and  that  no  plane  moves  unless  the 
loading  point  moves  it.  Under  such  an  assumption  the  yield  point 
in  compression  is  unaltered  by  loading  into  the  plastic  range  in 
tension.  This  however  is  not  in  accord  with  experiment.  Often 
hardening  in  tension  produces  a softening  in  compression  known 
as  the  Bauschinger  effect.  The  yield  points  in  tension  and  com- 
pression at  right  angles  to  the  original  direction  of  loading 
may  also  be  affected  (cross  effect).  In  theory  at  least  the 
whole  yield  surface  can  expand  uniformly.  In  general  of  coarse 
the  whole  yield  surface  is  affected  in  some  way  or  other  by  any 
hind  of  loading  and  not  just  part  of  the  yield  surface  as  was 
indicated  in  Fig.  1.  In  order  to  produce  these  effects  it  is 
necessary  to  assume  some  sort  of  interdependence  among  the  load- 
ing planes. 

The  number  associated  with  a plane,  which  completely 
determines  i**s  position,  is  the  distance  of  the  plane  from  the 
origin  (exactly  how  distance  is  to  be  defined  in  stress  space 
will  be  discussed  later).  It  seems  best  to  abandon  the  use  of  a 
loading  function  as  such  and  to  concentrate  attention  upon  the 
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changes  in  distance*  of  the  planes  from  the  origin  as  loading 
proceeds.  The  desired  interdependence  between  the  loading  planes 
may  be  achieved  by  postulating  a relation  between  the  motions  of 
the  various  planes.  It  will  be  assumed  that  the  loading  of  a 
given  plane  induces  a motion  in  any  other  plane  according  to 
some  rule  which  involves  the  orientations  of  the  two  planes  and 
the  measure  of  loading  of  the  given  plane.  Only  those  planes  in 
contact  with  the  loading  point  ray  be  said  to  be  loading,  thus 
the  motion  of  the  loading  point  completely  determines  the  total 
motion  of  the  planes  being  loaded.  However,  the  planes  being 
loaded  induce  a motion  in  each  other  which  when  subtracted  from 
their  total  motion  leaves  a remainder  to  be  accounted  for.  This 
remainder  which  is  called  the  direct  motion  will  be  taken  to  be 
the  measure  of  loading.  At  this  time  the  reader  is  asked  to 
accept  the  following  statement  on  faith.  If  a suitable  relation 
between  the  direct  and  induced  motions  is  given,  and  the  loading 
path  is  given,  then  the  requirement  that  the  direct  motion  shall 
always  be  in  the  loading  (outward)  sense  results  in  a determinate 
system.  A partial  justification  of  this  statement  is  made  later 
in  the  paper.  The  situation  is  a little  clearer  where  the  in- 
duced motion  is  very  small  for  then  the  yield  surface  is  almost 
the  same  as  though  the  Planes  moved  independently. 

There  is  at  least  an  intuitive  connection  between  these 
notions  and  the  behavior  of  crystal  aggregates.  During  plastic 


* 

Distance  in  stress  space  is  defined  later. 
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straining  the  inactive  slip  systems  are  hardened  as  well  as  the 
active  ones,  a phenomenon  known  as  latent  strain  hardening. 

This  suggests  the  induced  motion  of  loading  planes.  That  it  js 
possible  to  obtain  reasonable  results  basod  on  the  above  assump- 
tions has  been  shown  by  actually  constructing  a few  examples. 

In  these  examples  the  type  of  behavior  exhibited  by 
some  known  stress-strain  relations  possessing  a Dauschinger 
effect  [ 9 ] , [ 10]  , [ 11  ] has  been  duplicated.  The  yield  bound- 
ary after  biaxial  loading  according  to  a typical  theory  of  this 
kind  is  shov/n  in  Fig.  4.*  The  initial  yield  boundary  might  be 
the  J2  ellipse.  After  loading  in  tension  to  P this  ellipse  is 
enlarged  and  displaced  to  the  right.  In  an  example  constructed 
using  plane  loading  surfaces  acting  i.nterdependently  the  be- 
havior illustrated  in  Fig.  5 'was  obtained.  The  greater  part  of 
the  yield  boundary  is  modified  in  the  same  way  as  before  but 
now  a corner  appears  at  the  loading  point  P.  By  varying  certain 
parameters  involved  the  corner  at  P may  be  made  as  blunt  as  we 
please.  The  behavior  shov/n  in  Fig.  1 and  in  Fig.  4 may  be 
obtained  as  limiting  cases.  The  yield  boundary  according  to 
simple  J2  flow  theory  is  also  obtainable  as  a limiting  case. 

In  the  stress-strain  relations  for  this  theory  the 
direct  motion  of  a loading  plane  is  taken  as  the  measure  of  its 
contribution  to  the  plastic  strain  rather  than  the  total  motion 
so  that  for  a small  group  of  planes: 


Drucker's  [2]  definition  of  a work  hardening  material  does  not 
require  the  origin  to  bo  inside  the  yield  surface. 
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& e.^  = H(r-,,w)  n.  . 6u  (5) 

where  is  the  direct  motion.  The  integrated  form  of  (5) 
naturally  raises  the  question  of  path  independence.  About  the 
best  that  can  be  said  at  present  is  that  if  the  yield  surface 
determines  the  direct  motion  then  the  total  plastic  strain  is 
the  same  for  any  two  loading  paths  which  result  in  the  same  yield 
surface.  This  amounts  almost  to  a tautology.  The  difficulty 
is  that  although  the  direct  motion  determines  the  yield  surface 
the  converse  is  not  true  unconditionally.  Certainly  any  un- 
qualified assertion  concerning  path  independence  must  be  false 
in  some  of  the  limiting  cases.  However  there  are  cases  in  which 
path  independence  is  possible. 

Figure  6 illustrates  a case  based  on  the  examples 
constructed  in  a later  section  in  which  two  different  loading 
paths  OPQ  and  0P*Q  result  in  the  same  yield  surface  and  total 
plastic  strain.  There  is  a region  of  total  loading  beyond  Q 
into  which  a nearly  radial  loading  path  may  proceed;  it  is  shown 
bounded  by  the  dashed  lines.  In  the  limiting  cases  in  which  the 
corner  becomes  blunted  the  region  of  total  loading  narrows  down 
to  a line  so  that  nearly  radial  loading  becomes  strictly  radial 
loading. 

The  mathematical  formulation  of  the  problem  of  finding 
the  yield  surface  according  to  the  present  theory  involves  an 
integral  equation.  The  conditions  under  which  the  problem  has 
a unique  solution  have  not  yet  boon  found.  There  is  no  doubt 
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about  uniqueness  in  the  example  solutions  given  but  only  nearly 
radial  loading  paths  have  been  considered.  Even  if  suitable 
conditions  are  imposed  to  insure  the  unique  determination  of 
the  direct  motion  it  remains  to  be  shown  that  the  corresponding 
stress-strain  relations  lead  to  unique  solutions  of  boundary 
value  problems.  In  order  to  extend  the  uniqueness  theorem  to 
cover  the  case  of  interdependent  motion  it  may  be  that  some 
additional  restrictions  must  be  placed  on  the  rule  governing  the 
relation  between  the  direct  motion  and  the  induced  motion. 

The  problem  of  determining  the  vield  surface  for 
loading  paths  which  are  not  nearly  radial  is  considerably  more 
complicated  than  in  the  case  in  which  the  loading  planes  move 
independently.  The  yield  surface  resulting  from  such  a loading 
path  might  look  something  like  Fig.  7 in  which  OPQ  is  a broken 
line  with  a sufficiently  large  turning  angle  at  P.  A corner 
previously  formed  at  P has  boon  rounded  off  and  a new  corner 
formed  at  Q. 

The  present  theory  affords  a means  of  taking 
Bauschinger  and  allied  effects  into  account  without  explicitly 
introducing  strain  into  tho  loading  functions  as  has  been  the 
usual  practice  to  date.  However  tho  attendent  complications  in 
determining  the  yield  surface  have  not  boon  avoided.  Although 
It  is  still  convenient  to  spoa?'  of  loading  pianos  there  doesn't 
seer-!  to  be  anything  left  which  could  strictly  be  called  a loading 
function.  The  direct  motion  is  sort  of  an  intermediate  variable 
determined  by  the  loading  history  and  which  in  turn  determines 


Foitcr's  theorem  [1]  is  applicable  in  the  case  of  independent 
motion. 
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both  the  yield  surface  and  the  plastic  strain. 

Stress-Strain  Relation  for  Isotropic  Materials . 

In  this  section  the  loading  planes  are  assumed  to 
move  independently.  A general  stress-strain  relation  for  iso- 
tropic materials  is  constructed.  The  stress-strain  relations 
for  slip  theory  and  for  a theory  which  agrees  v/ith  J2  deformation 
theory  for  nearly  radial  loading  are  written  down  as  special  cases. 

In  dealing  with  plane  loading  surfaces  a certain 
amount  of  mathematical  machinery  of  a geometrical  nature  is 
almost  indispensable  for  concise  expression  of  the  ideas  involved. 
Once  a suitable  stress  space  has  been  constructed,  and  some 
appropriate  coordinate  system  has  been  defined  in  this  space, 
writing  down  a stress-strain  relation  for  isotropic  materials 
based  on  infinitely  many  plane  loading  surfaces  is  almost  a 
trivial  matter.  Since  plastic  strain  is  assumed  to  depend  only 
on  the  stress  deviator,  which  has  five  independent  components, 
a stress  space  of  five  dimensions  is  indicated.  For  the  five 
cartesian  coordinate  members  it  is  possible  to  choose  linear 
combinations  of  the  stress  deviator  such  that  a ’’stress  vector" 
in  stress  space  transforms  liVe  a vector  when  the  corresponding 
stress  tensor  is  transformed  by  a rotation  of  the  physical  co- 
ordinates. The  length  of  the  stress  vector,  being  an  invariant, 
is  not  changed  by  the  transformation.  This  inevitably  leads  to 
a distance  function  in  stress  space  proportional  to  (Jp)  • In 
such  a stress  space  any  quantity  such  as  distance,  angle,  area, 

* J2  = is±j  stj  . 
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etc.  defined  in  terms  of  the  metric  may  be  used  consistently 
because  it  is  an  invariant  v/ith  respect  to  chances  of  coordinates 
in  physical  space.  A stress  space  of  this  kind  is  constructed 
in  a’ipe nd  ix  A.  Several  systems  of  coordinates  for  this  stress 
space  are  defined  in  the  appendix  which  also  contains  other  de- 
tails of  a geometrical  nature.  Only  a brief  description  of  these 
results  will  be  given  in  this  section  as  the  need  arises. 

The  contribution  to  the  total  plastic  strain  of  a 
group  of  planes  which  were  initially  tangent  to  the  initial  yield 
surface  over  an  element  of  area  da  is  assumed  to  be  given  by: 

6 ea  = H(f,na)  na  da  (a  = 1,2,.., 5)  (6) 

where  ea  is  the  plastic  strain  vector,  and  na  is  the  unit  normal 
vector  to  the  element  of  area  directed  in  the  loading  sense.  The 
stress  vector  is  denoted  by  xa  and  f is  given  by  f = xQ  , 
or  rather  the  maximum  value  of  this  quantity  attained  during  the 
course  of  loading.  H > 0 if  f > Ra  where  Ra  is  the  vector 

from  the  origin  to  a point  within  the  element  of  area  da  on  the 
initial  yield  surface,  otherwise  H vanishes.  The  contributions 
of  all  planes  are  summed  to  get  the  total  plastic  strain: 

ea  = J II  na  da  (7) 

S is  the  initial  yield  surface  v/hich  is  the  envelope  of  the 
initial  positions  of  all  the  loading  planes.  The  stress-strain 
relation  (7)  is  in  general  anisotropic;  further  specialization 
is  necessary  to  arrive  at  a relation  applicable  to  initially 
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isotropic  materials. 

In  appendix  A a special  spherical  coordinate  system 
is  defined  in  which 

J2  ' P2  = xa  xa 

and  (8) 

J = JL 

3 ^.n/3 

Three  other  coordinate  angles  6^,  ©2,  ©3  are  introduced  and  the 
transformation  of  coordinates  is  given  by  equations  of  the  form: 

xa  = p "a  (V  (P  = M>3>4>  (9) 

where  KaKa  = 1.  The  equation  of  the  initial  yield  surface  of 
an  isotropic  material  in  these  coordinates  is  given  by: 


cos 


3©4 


p=  p<v 


(10) 


Isotropy  requires  that  all  planes  be  treated  equally 
which  are  initially  tangent  to  the  initial  yield  surface  along 
its  intersection  with  a coordinate  surface  ©4  = constant.  The 
special  form  of  the  stress-strain  relation  (7)  for  initially 
isotropic  materials  is  thus: 

0 


F ( f , ©4 ) na  da 
(js 


(11) 


or  more  explicitly: 
0 


ea  = 


- |r  sin  e2  sin  3e4  dex  de2  d®3  de4  • 


(12) 
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where  H = 0,  f < P na  £a  . Sometimes  some  other  coordinate 
system  is  more  convenient  than  the  one  used  in  (11);  in  a 
general  system  of  coordinates  this  relation  may  be  written: 

ea  = j H(f,J3*)  na  da  (13) 

where  J-^*  is  the  value  of  at  the  point  of  tangency  to  the 
yield  surface  of  the  loading  plane  in  its  initial  position. 

The  simplest  special  case  of  the  stress-strain  relation 
constructed  above  is  that  in  which  the  initial  yield  surface  is 
p = k and  all  planes  are  treated  equally,  in  this  case; 

■i 

ea  = H(f)  <a  da  (H  = 0,  f < k)  (l4) 

els 

The  initial  yield  surface  is  of  course  the  Mises  yield  surface 
which  is  a sphere  in  the  present  stress  space.  For  radial 
loading  (or  nearly  radial  loadin'-  for  that  matter)  the  strain 
vector  must  have  the  same  direction  as  the  stress  vector  because 
of  the  perfect  symmetry  of  the  conf iguration  of  loading  planes. 
Again  because  of  symmetry  the  length  of  the  strain  vector  can 
only  depend  on  the  length  of  the  stress  vector  and  not  upon  its 
direction.  The  stress-strain  relation  for  nearly  radial  loading 
must  be: 


ea  = Ftx^xp)  xa  05) 

In  the  more  usual  notation  this  is 

Cij  = F(J2)  sij  (l6> 
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which  is  immediately  recognized  ?s  the  stress-strain  relation 
of  Jo  deformation  theory.  The  form  of  the  stress-strain  rela- 
tion for  a broken  line  loading  path  with  a large  turning  angle 
is  given  in  appendix  B. 

Another  interesting  special  case  is  slip  theory.  The 
loading  planes  of  slip  theory  in  their  initial  positions  are 


given  byt 


f = = h sjL.  (\t  ^ = k 


(17) 


where  and  are  all  possible  pairs  of  orthogonal  unit  vectors 
in  physical  space  [l  ] . It  is  easy  to  verify  that  the  vector t 


sij  ' ,c(h  'p  * h -“d 


(18) 


drawn  from  the  origin,  is  normal  to  the  corresponding  plane 

(in  9 dimensional  stress  space)  and  moreover  terminates  in  the 

plane.  Also  £ s^j  sij  = k^  and  ^ s^?  s^°  s^°  = 0.  Therefore 

the  planes  of  slip  theory  are  all  those  planes  tangent  to  the 
2 

surface  J2  = k'  where  = 0 and  no  others.  The  same  holds 
in  x stress  space.  All  these  planes  are  to  be  treated  equally 
because  slip  theory  applies  to  an  initially  isotropic  material. 
Here  is  a case  in  which  not  all  planes  available  in  stress  space 
are  put  to  use.  Even  though  all  the  planes  used  are  initially 
tangent  to  the  sphere  p = k the  initial  yield  surface,  as  is 
well  known,  is  the  Tresca  yield  surface.  The  stress-strain 
relation  is  obtained  from  (12)  by  omitting  the  integration  with 
respect  to  9^. 
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ea  = | H(f)  sin  G2  A01  d©2  d©3  (H  = 0,f  < k)  (19) 

where  S'  is  the  intersection  of  J2  = k with  = 0.  The 
parameters  9i»®2j®3  involved  in  the  original  formulation  of 
slip  theory  [5]v/ere  somewhat  different. 

In  both  of  the  preceding  cases  the  function  H(f)  may 
be  obtained  from  a stress-strain  curve  by  solving  (in  closed 
form)  a simple  type  of  integral  equation.  In  the  general  case 
more  experimental  information  than  a stress-strain  curve  would 
be  needed  to  determine  H^Q^),  and  solving  the  corresponding 
integral  equation  would  not  be  a simple  matter.  Of  course  a 
form  for  H could  be  assumed  involving  a number  of  arbitrary 
parameters  and  then  one  could  attempt  to  match  the  data  by 
adjusting  the  parameters,  but  at  present  it  doesn't  seem  to  be 
worth  while  to  pursue  the  matter  any  further. 

If  the  loading  path  is  such  that  the  principal  axes 
of  stress  do  not  rotate  then  the  loading  path  is  confined  to 
the  x-p  x2  plane  provided  the  physical  coordinates  are  referred 
to  the  principal  axes.  Again  if  the  only  stresses  imposed  are 
say  d-^i  and  cs^2  then  the  loading  path  is  confined  to  the  x-^x^ 
plane.  In  these  cases  the  stress-strain  relations  based  on 
loading  planes  in  stress  space  can  be  reduced  to  a form  having 
the  appearance  of  a stress-strain  relation  based  on  "loading 
lines"  in  a "stress  plane".  The  reduction  can  always  be  made, 
at  least  in  theory.  For  simplicity  the  argument  will  be  pre- 
sented only  for  the  case  in  which  the  initial  yield  surface  is 
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p = k and  for  the  sake  of  being  definite  suppose  the  loading 
path  is  confined  to  the  x-l,x2  plane. 

Consider  a plane  which  has  the  trace  t in  the  x-^x^ 
plane  before  loading  begins  (see  Fig.  8).  Let  OP  be  the 
perpendicular  from  0 to  t and  let  9^  be  the  angle  between  OP 
and  the  x-^  axis.  Let  OQ  be  the  perpendicular  from  0 to  the 
plane  in  question  and  let  ©2  be  the  angle  POQ,  Suppose  that 
during  the  course  of  loading  this  plane  is  pushed  out  so  that 
its  trace  is  now  t1 . The  distance  of  the  plane  from  the  origin 
which  was  OQ  = k is  now  OQ1  = OP7  cos  ©2  = Consider  now  a 
spherical  coordinate  system  in  which  the  x^  and  x2  coordinates 
of  Q are  given  by: 


Xt  = k cos  0o  cos  9, 


x2  = k cos  92  sin  9-^ 


(20) 


and  x^jX^x^  are  given  in  terms  of  9-p©2  and  two  more  parameters 
(see  Per  example  appendix  A Eq.  A-28).  Now  since  there 
is  a two  parameter  set  of  points  Q on  the  sphere,  all  with  the 
same  9-^  and  ©2  coordinates,  there  is  a two  parameter  family  of 
planes  tangent  to  the  sphere  with  the  same  trace  t in  the  x-^,x2 
plane  as  the  plane  considered  above.  The  contribution  of  these 
planes  to  the  plastic  strain  is  given  by: 


62eQ  = 


2k  /iTt 


=■ 


p = J H(f,J3  ) sin-  292  sin  9^  d9^  d9^  d9^  d©2 


= H (f,e1,©2)  ^ sin-  2©^  d©1  d©2  (p  r,  1,2)  (21) 


All-106 


19 


where  Xj_  = cos  ©-^  , £2  = s^n  ®1  an<^  H = 0,  f < k.  For  values 
of  ©2  between  0 and  x/2  there  is  a continuous  distribution  of 
parallel  traces  each  representing  a two  parameter  family  of 
planes. * To  obtain  the  strain  due  to  moving  out  all  those 
planes  whose  traces  fall  between  the  one  tangent  to  the  circle 
and  the  one  at  t'  requires  the  integration  of  6*TSp  with  respect 
to  9p  between  the  limits  0 and  cos”-*-  (k/OP'  ).  Lot  OP*  = f,  then 

Ucos"1(k/f)  _ P 

6 ep  = j = j H(f  cos  02)6^)  ©2  ) K p sin  ^ ©2  d ©2  d ©-^ 

= I (f,01)  ^ d©1  ; I = 0,  f < k ; ( p = 1,2)  (22) 

6 ep  may  be  interpreted  as  the  strain  due  to  moving  out  a 
"loading  line"  initially  tangent  to  the  yield  boundary,  moreover 
6ep  is  normal  to  this  loading  line.  So  long  as  the  loading  is 
in  the  x^,x2  piano  6e^  = 6e^  = 6e^  = 0 for  an  isotropic  material, 
otherwise  the  remaining  strains  rnay  be  expressed  as  single 
integrals  if  so  desired.  Note  that  up  to  this  point  the  exact 
shape  of  the  (plane)  loading  path  has  not  entered  into  the 
calculations.  To  complete  the  analogy  the  plastic  strain  ep 
is  given  as  the  sum  of  the  contributions  from  all  the  loading 
lines . 

02ri-  _ 

e0  = \ H (f,©1)  d©1  * K = 0,  f < k ; (p  = 1,2)  (23) 

■ 0 o 

In  appendix  B the  foregoing  reduction  is  carred  out  for  the 
case  in  which  all  planes  arc  treated  equally.  A similar  reduc- 
tion to  two  dimensional  loading  planes  in  a three  dimensional 
stress  space  could  be  made. 

T 


The  trace  tangent  to  the  circle  represents  one  plane  only. 
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Induced  Motion, 

In  the  first  section  it  was  pointed  out  that  a 
realistic  theory  of  stress-strain  relations  based  on'  plane 
loading  surfaces  should  provide  for  some  sort  of  interdependent 
motion  of  planes  to  properly  account  for  the  changes  in  shape 
of  the  yield  surface  as  loading  proceeds.  A few  assumptions 
regarding  the  nature  of  this  interdependence  were  made.  Loading 
of  any  given  plane  was  assumed  to  cause  an  induced  motion  in  all 
other  planes.  It  was  proposed  to  separate  the  motion  of  a plane 
into  two  parts,  namely  direct  and  induced,  where  the  direct 
motion  is  associated  with  the  loading  of  a given  plane  and  the 
induced  motion  represents  the  effects  on  the  given  plane  of 
loading  on  all  other  planes.  In  this  section  enough  analytical 
detail  is  supplied  to  apply  the  theory  to  a few  simple  cases. 

The  contribution  to  the  induced  motion  of  a given 
plane  due  to  the  direct  motion  of  an  arbitrary  small  group  of 
planes  is  assumed  to  be: 


6 r2(0)  = F [ r1(6'  ),  6,  9'  ] da  (r1  > 0) 


(24) 


=0  (r1  = 0) 

where  Q and  O'  denote  the  orientations  of  the  given  and  arbitrary 
planes  respectively.  The  total  induced  motion  of  the  given 
plane  is  obtained  by  summing  the  contributions  from  all  planes: 


r„  = 


F(r1,  9,  Q»)  da 


(25) 
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The  total  distance  of  any  loading  plane  from  the  origin  is 
given  by: 

r = rQ  + ri  + r2  (26) 

where  r is  the  initial  distance  from  the  origin.  In  the  first 
o 

section  certain  other  assumptions  were  made  concerning  the 
direct  motion.  These  were  that  only  those  planes  moving  in 
contact  with  the  loading  point  could  receive  any  direct  motion 
and  that  direct  motion  is  a non-decreasing  function  of  time  as 
loading  proceeds.  In  symbols  these  are: 


r - n<x 

xa  > 

0 

=$>  rx  = 0 

(27) 

• 

rl 

> 

0 

(28) 

It  is  understood 

that: 

r 

- % 

Xa 

> 0 

(29) 

Note  that: 

r - n<x 

• 

*d  = 

0 

^ > 0 

(30) 

In  addition  to  these  assumptions  it  is  expected  that  certain 
restrictions  must  be  placed  upon  F in  order  to  obtain  a con- 
sistent theory. 

Two  examples  are  constructed  in  the  following  to 
check  the  general  reasonableness  of  the  assumptions  so  far  and 
to  gain  some  insight  into  the  nature  of  the  problem  of  deter- 
mining the  yield  surface  according  to  the  present  scheme.  To 
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keep  the  calculations  as  simple  as  possible  the  initial  yield 
surface  is  taven  to  be  p = 1 and  the  form  of  (25)  is  taken  to 
be : 


= 


r^(X  + p cos  iO  da 


(3D 


where  X and  |i  are  constants  and  4)  is  the  angle  between  the 
normals  to  the  planes  whose  orientations  were  denoted  by  0 and 


6'  . 

In  the  first  example  take  |i  = 0,  X > 0 and  let  the 
loading  path  be  coincident  with  the  positive  x-^  axis.  From 
(31)  r2  is  a non-decreasing  function  of  time  alone.  Thus  the 
yield  surface  expands  uniformly  except  for  the  portion  of  the 
surface  affected  by  direct  motion.  As  soon  as  the  loading  point 
reaches  x-^  = 1 some  planes  begin  to  lead.  Assume  for  the  moment 
that  once  a plane  begins  to  load  it  does  not  subsequently  unload. 
Then  at  any  instant  all  Planes  which  have  been  loaded  pass  through 

the  loading  point  and  all  planes  which  have  not  been  loaded  are 
tangent  to  the  expanding  spherical  part  of  the  yield  surface; 
the  planes  just  beginning  to  load  do  both.  The  picture  is  as 
in  Fig.  9.  A corner  in  the  yield  surface  appears  at  the  loading 
point  P.  The  radius  of  the  spherical  part  of  the  yield  surface 
as  a function  of  the  load  is  calculated  in  appendix  0.  As 
X oo  the  corner  at  P becomes  blunter  and  vanishes  in  the  limit. 


In  the  appendix  it  is  shown  that  for  X finite  the  loading  path  may 
turn  through  a certain  angle  without  unloading  any  planes  which 
have  already  been  loaded.  Suppose  the  loading  point  follows 


7 
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such  a bent  path,  then  at  any  instant  the  situation  is  exactly 
the  same  as  before  except  that  P lies  a little  ways  off  the  x-^ 
axis.  The  same  configuration  of  loading  planes  would  have  been 
produced  had  the  loading  point  reached  P by  a radial  path  so 
the  direct  motion  of  any  nlane  is  the  same  in  either  case.  As 
X —)  00  the  allowable  turning  angle  goes  to  zero  so  there  is  no 
path  independence  in  that  case.  If  X = 0 then  there  is  no  in- 
duced motion  and  we  have  the  case  of  the  last  section.  If  X 
is  negative  there  is  the  possibility  of  two  configurations  of 
loading  planes  for  a given  value  of  the  load.  Proper  restric- 
tions on  F should  rule  out  negative  values  of  X , at  least  when 

p = 0. 

For  the  sjeond  example  take  X = 0 and  p > 0;  (31)  be- 
comes 

n p 

T2  = p 1 r-j_  cos  \jj  da  = p r,  na  na'  da 

(Is  Js 

= na  | rl  V da  = na  Aa  (32) 

where  Aa  are  constants.  There  is  a simple  geometrical  inter- 
pretation of  the  induced  motion  in  this  case.  Suppose  the 
yield  surface  has  attained  some  irregular  shape  through  loading 
and  then  additional  loading  takes  nlace.  What  is  the  change  in 
shape  of  the  yield  surface  due  to  the  additional  induced  motion 
alone?  Consider  a plane  tangent  to  the  ”-ield  surface  at  some 
arbitrary  point  P.  Suppose  P is  displaced  to  F'  by  the  displace- 
ment vector  A a carrying  the  plane  along  with  it  (see  Fig.  10). 
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The  distance  the  plane  has  been  displaced  parallel  to  itself  is 
A<j  . Thus  the  induced  motion  given  by  (32)  has  the  effect 
of  displacing  the  whole  yield  surface  like  a rigid  body.  For 
pL  positive  the  displacement  is  in  the  same  general  direction  as 
the  additional  leading.  Of  course  the  direct  motion  further 
distorts  the  yield  surface  in  the  vicinity  of  the  loading  point. 
The  behavior  for  radial  loading  is  similar  to  that  of  the  first 
example,  that  is  for  |j.  finite  a corner  appears  at  the  loading 
point  and  again  nearly  radial  loading  paths  exist.  In  this 
case  also  the  corner  disappears  for  |i  = co  and  nearly  radial 
becomes  strictly  radial.  The  calculations  are  given  in  appendix 
C.  For  sufficiently  large  negative  values  of  y.  there  is  the 
possibility  of  two  configurations  of  loading  planes  for  a given 
value  of  the  load.  The  reader  may  notice  other  inconsistencies 
for  negative  values  of  X and  jj.  ; at  any  rate  these  examples  show 
that  there  is  a definite  need  for  restrictions  on  F, 

A combination  of  these  two  examples  gives  the  type  of 
behavior  referred  to  in  the  first  section.  Here  again  the  stres 
strain  relation  would  agree  with  Jo  deformation  theory  for 
nearly  radial  loading  paths. 

Observations  and  Conclusions . 

The  plane  loading  surface  has  been  used  as  the  funda- 
mental building  block,  so  to  speak,  in  the  construction  of  stres 
strain  relations  of  a quite  general  type.  Considerable  flexi- 
bility has  been  gained  by  assuming  the  displacements  of  the 
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various  planes  to  be  interrelated.  In  many  special  cases  of  the 
present  theory  the  plastic  strains  are  path  independent  to  a 
certain  extent  and  a special  class  of  loading  paths  exist,  called 
nearly  radial,  for  vhich  the  stress-strain  relations  reduce  to 
those  of  A deformation  theory.  Nearly  radial  loading  paths  were 
shown  to  exist  in  all  cases  not  involving  induced  motion  and 
even  in  some  cases  where  induced  motion  of  the  loading  planes  is 
allowed.  It  is  quite  probable  that  nearly  radial  loading  paths 
exist  whenever  the  function  F is  bounded.  Further  theoretical 
investigations  should  be  made  to  settle  the  question.  Certainly 
the  validity  of  the  theory  should  be  tested  by  comparison  with 
experiment. 

It  was  remarked  in  an  earlier  section  that  the 
yield  surfaces  of  some  known  theories  could  be  duplicated 
by  limiting  cases  of  the  present  theory.  It  is  extremely 
interesting  to  look  into  the  corresponding  question  con- 
cerning the  two  stress-strain  relations.  Only  a cursory 
examination  is  necessary  to  show  that  the  known  stress  - 
strain  relations  cannot  bo  obtained  as  limiting  cases  of 
the  present  theory.  Consider  the  case  of  the  first  example 
discussed  in  the  text  and  let  X be  very  large;  then  the 
yield  surface  is  very  nearly  that  according  to  simple  J2 
flow  theory.  Suppose  loading  in  tension  has  proceeded  up 
to  a certain  point  and  then  the  load  is  removed.  Now 
compare  the  stress-strain  curves  for  loading  in  com- 
pression according  to  J2  flow  theory  and  according  to  the 
present  theory.  According  to  J2  flow  theory  there  will  be 
a sharp  break  in  the  stress-strain  curve  at  the  yield  point, 
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but  according  to  the  present  theory  the  stress-strain  curve  will 
be  smooth  because  the  planes  tangent  to  the  sphere  on  the  com- 
pression side  have  not  been  previously  loaded.  At  first  sight 
this  may  seem  distressing,  one  might  have  hoped  that  the  present 
theory  would  be  general  enough  to  include  sue’1,  a classic  example 
as  Jg  flow  theory  as  a limiting  case.  However,  the  experimental 
facts  are  in  favor  of  the  present  theory,  qualitatively  at  least. 
Stress-strain  curves  in  compression  following  stressing  in  tension 
are  always  found  to  be  quite  smooth  [ 12 ] . Even  those  linear 
incremental  stress-strain  relations  which  exhibit  a Bauschingcr 
effect  predict  a sharply  defined  yield  point  for  such  an  experi- 
ment as  that  described  above. 

More  complicated  experiments  have  been  performed  in 
which  the  material  is  first  stressed  in  tension  and  then  loaded 
in  some  way  other  than  mere  compression.  The  experiments  of 
Taylor  and  Quinney  [13  ] and  of  Klingler  and  Sachs  [ 14]  were  of 
this  type.  In  all  cases  the  reported  stress-strain  curves  for 
the  second  loading  are  quite  smooth  as  would  be  predicted  by  a 
stress-strain  relation  of  the  type  considered  in  this  paper.  On 
the  other  hand  it  is  a well  known  fact  that  the  stress-strain 
curve  for  a second  loading  in  tension  does  have  a rather  sharply 
defined  yield  point.  This  singular  behavior  at  a point  on  the 
yield  surface  where  loading  has  recently  occured  may  be  regarded 
as  evidence  for  corners  in  the  yield  surface.  The  fact  that  any 
intermediate  loading  (before  the  second  tension  test)  is  likely 
to  destroy  the  sharpness  of  the  yield  point  is  also  evidence  in 
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support  of  a theory  of  the  present  type.  For  some  reason 
(possibly  creep)  the  stress-strain  curve  for  reloading  in  tension 
does  not  always  have  a perfectly  sharp  break  at  the  yield  point; 
the  same  reason  may  explain  why  corners  have  not  always  been 
found  by  those  experimenters  who  have  looked  for  them  [ 8 ] , [ 15  ]. 

Certainly  the  present  theory  is  far  from  being  well 
developed  and  much  remains  to  bo  done  before  it  can  even  claim  to 
be  acceptable.  Sven  so  fundamental  a condition  as  the  work 
hardening  condition  is  not  automatically  satisfied  unless  re- 
strictions arc  placed  on  the  rule  governing  the  induced  motion. 
These  restrictions,  and  possibly  others  necessary  to  insure 
uniqueness,  are  yet  to  be  discovered. 
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Appendix  A 


Five-dimensional  stress  space . 

If  plastic  strains  are  assumed  to  be  incompressible 
then  the  stress-strain  relations  and  the  yield  surface  depend 
only  on  the  stress  deviator  and  not  upon  = dj^  . This  means 
that  in  the  usual  nine-dimensional  stress  space  the  yield  sur- 
face is  a cylinder  normal  to  the  plane  dj^  = 0.  Rather  than 
represent  stress  in  this  nine-dimensional  space  in  which  each 
component  of  stress  acts  as  a coordinate  it  is  more  convenient 
for  present  purposes  to  choose  some  other  system  of  coordinates 
which  better  fits  the  yield  surface.  Since  the  stress  deviator 
has  only  five  independent  components  it  is  possible  to  represent 
it  geometrically  in  a suitably  chosen  space  of  five  dimensions. 
Methods  similar  to  the  one  given  below  for  constructing  a stress 
space  exist  in  the  literature  [l6  ] , [17  ] but  they  are  usually 
not  given  in  much  detail  so  the  one  used  in  this  paper  will  be 
discussed  in  this  appendix. 

The  stress  space  chosen  here  is  a subspace  of  the 
usual  nine-dimensional  stress  space.  The  coordinates  in  this 
subspace  are  chosen  in  the  following  way.  The  x^  axis  (which  will 
later  be  ignored)  is  taken  normal  to  the  plane  dj_j_  = 0.  The 
axis  is  the  projection  of  the  o-q  axis  on  the  plane  d^  = 0 
and  the  Xp  axis  is  taken  normal  to  the  x-j_  and  x^  axes  so  that 
x^,  Xp,  x^  forms  a right-handed  system  of  coordinates.  The 
scale  of  measurement  on  the  x-^,  Xp , x^  axes  is  distorted  by  a 
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factor  of  for  a reason  which  will  be  given  later.  The  equa- 
tions of  transformation  between  the  two  systems  of  coordinate  are; 


X1  = 
x2 


— — ■ (2d-i-t  •*  Cps  *■  dpp) 

2/3  33 


= I (d22  - d^) 


( A-l) 


x^  = 


— ^ - (d^-^  °22  d3^ ^ 


Three  more  axes  x^,  x^,  and  x^  are  added  to  accommodate  shears: 


X3  ~ T12 


X)+  = t 


23 


(A-2) 


x5  = T 13 


so  that  x^,  x2  x^  is  a rectangular  cartesian  system  cf  co- 

ordinates. The  reason  for  the  factor/?  is  the  following:  under 
a rotation  of  physical  coordinates  d jj  behaves  like  a second 
order  cartesian  tensor,  the  factor  f?  is  necessary  in  order  that 
xa  , (a  = 1,2,... 6),  should  transform  like  a cartesian  vector. 

The  relations  (1)  and  (2)  may  be  compactly  written  as 

follows : 


xa  - dij 


(A-3) 


in  which  3a,  ^ are  a set  of  3x3  matrices  whose  definitions  are: 


s 


I 
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Since  all  formulas  concerning  plastic  strains  do  not  involve  the 
coordinate  this  coordinate  will  be  dropped  and  the  range  of 
summation  on  the  Greek  indices  will  hereafter  be  understood  to 
run  from  1 to  5.  The  quantities  satisfy  the  following 

identities  which  are  easily  verified. 


Ba’  ij  • Ba’  ji 
Ba’ij  = ^ 6ap 


B 


a>  ii 


= 0 


B 


a’lj 


B 


a’k't 


1 

5 


<& 


ik  V + 6i^6jk  ~ | 6ij  &k£  }* 


The  components  of  xa  depend  only  on  the  stress  devistor  as  is 
shown  directly  by  the  following: 


x a _ B d,  . 6j  . - B 


10  ij  cc’ij  10 


-i.1  (sij  + s V - Vlj  sij 


+ s Ba’ii  Ba’ij  Si 


ij 


(A  -6) 
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Equation  6 may  be  solved  for  the  components  of  the  deviator  in 


1 


/ 


! * 
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= 2 Boc>ij  xa 


( A-7) 


From  (7)  and  (5)  it  follows; 


J2  = * s13s13  = 2 xaxp  = = xaxa 


(A-8) 


that  is  the  length  of  a stress  vector  xa  is  2* 

Analogous  to  the  stress  vector  a plastic  strain  vector 
is  defined  as  follows: 


= 2 Ba)i;j  e.J  (A-9) 

The  stress-strain  relation  for  a linear  stress  theory  of  plas- 
ticity reads  as  follows  in  the  present  notation: 


de_  = G(f) 


a 


3f 


df 


(A-10) 


as  may  easily  be  verified.  The  tensor  character  of  this  relation 
is  preserved  in  the  present  representation  and  as  before  the 
plastic  strain  increment  is  normal  to  the  loading  surface. 

A coordinate  system  for  isotronic  materials . 

The  geometrical  representation  of  J?  as  given  in  the 
last  section  is  quite  simple  but  what  about  J^?  Explicitly  is: 


= is.,  . s .n.Si.-t  = x-,3  ^ 


'3  = rij5^5ki 


3 & 


JL 

\fz 


Xl(x22  + Xh2)  + A Xl(x32 
+ x^2)  + x^Cx  2 - x^2)  + 2x-,xLxc. 


(A-ll) 


\ 

*5 
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which  is  a rather  awkward  expression.  There  is  a simple  repre- 
sentation of  in  case  Sj_j  = 0,  i 4 j,  in  which  case  x^  = x^ 

= x^  = 0,  If  polar  coordinates  are  introduced  into  the  1,2  plane 
letting : 


x^  = p cos  9 
X2  = p sin  9 


then  we  find: 


= (P  ; J„  = 2j^ 


cos  39 


3tf  3 


The  corresponding  values  of  s . . are: 


S11 


= t=£-  cos  9 

/3 


s22  = - cos  (©  + |) 

s = - iR  COS  (9  - 

33  \/3  i 


(A-12) 


(A-13) 


(A-l4) 


The  reason  for  introducing  J2  and  J-  is  of  course  that  they  are 
invariant  under  rotations  of  the  physical  coordinates  and  hence 
play  a prominent  role  in  the  description  of  initially  isotropic 
materials.  V.'e  would  like  to  construct  a coordinate  system  in 
stress  space  convenient  for  such  a description.  In  order  to  gain 
a little  insight  into  the  problem  we  consider  the  effect  on  the 
vector  xa  of  subjecting  the  corresponding  tensor  Sjj  to  a rotation 
in  physical  space.  As  is  easily  checked  using  (5)  the  new  vector 
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xa'  corresponding  to  the  new  tensor  s^j  is  related  to  the  vector 
xa  by  a rigid  body  rotation  in  stress  space.  That  is  to  say 


V = A 


a?  XP 


(A-15) 


where  Aap  is  an  orthogonal  matrix.  If  is  the  rotation 

tensor  in  physical  space  then  Aap  is  given  by: 


la[3 


= 2 a 


ik  ai£  Ba’ij  BP’ k^t 


(A-l6) 


Thus  to  every  rotation  in  physical  space  there  corresponds  a 
rotation  in  stress  space.  The  correspondence  between  rotations 
however  is  not  one  to  one.  There  is  a ten  parameter  group  of 
rotations  in  the  five  dimensional  stress  space  but  only  a three 
parameter  sub-group  of  them  corresponds  to  the  three  parameter 
group  of  rotations  in  physical  space.  The  somewhat  messy 
expression  (11)  is  of  course  an  invariant  of  this  sub-group. 

The  correspondence  between  rotations  in  physical  space 
and  a certain  group  of  rotations  in  stress  space  suggests  at 
least  one  way  in  which  a coordinate  system  especially  adapted 
to  the  description  of  isotropic  materials  may  be  defined.  We 
first  note  that  for  a given  stress  tensor  the  coordinate  system 
in  physical  space  can  always  be  rotated  into  a system  of  princi- 
pal axes  for  that  stress  tensor.  This  means  that  in  stress 
space  an  arbitrary  vector  xa  can  always  be  rotated  into  the  1,2 
plane  by  a rotation  belonging  to  the  above  mentioned  sub-group. 
This  can  be  accomplished  in  six  different  ways  corresponding  to 
the  six  choices  of  right  handed  principal  axes  in  physical  space, 
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This  means  that  in  stress  space  an  arbitrary  vector  xQ  can  be 
rotated  into  any  one  of  six  equal  sectors  in  the  1,2  plane  by  a 
rotation  belonging  to  our  sub-group.  Thus  the  stress  vectors 
xn  in  the  sector  of  the  1,2  plane  0 < © < 2.  are  representative 
of  all  vectors  in  stress  space  through  these  rotations.  The 
coordinate  system  is  constructed  as  follows.  Choose  a vector 
in  the  1,2  plane  given  by  its  length  p and  the  angle  ©^  it  makes 
with  the  x^  axis  where  0 < ©^  < 2.  . Next  apply  an  arbitrary 
rotation  belonging  to  our  sub-group  which  will  involve  three 
parameters,  call  them  6-p  ©2  and  0^.  The  resulting  vector: 


xa  ~ xa  ®2  ’ ®3  ’ 


(A-17) 


can  be  any  vector  in  stress  space  so  (16)  can  be  regarded  as  the 
equations  defining  a transformation  from  x coordinates  to  p , © 
coordinates.  With  suitable  restrictions  on  the  range  of  S^,  ©2 
and  0-^  the  new  coordinate  system  covers  stress  space  once  and 
only  once.  Obviously  the  values  of  J2  and  are  the  same  for 
any  two  stress  states  which  have  the  same  p and  ©^  coordinates 
because  these  quantities  are  unaffected  by  rotations.  Expres- 
sions for  these  invariants  are  still  given  by  (13)  where  we  take 
© to  be  ©^.  The  principal  deviatoric  stresses  in  the  order 
S1  2 sp  ^ s3  aro  given  by  (14)  when  ©^  is  restricted  to  the 
range  0 ©^  < 2 . As  an  example  of  the  use  of  this  coordinate 

system  note  that  the  equation  of  the  initial  yield  surface  of 
an  isotropic  material  must  bo  of  the  form: 
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p = p (el(.)  (a-i8) 

The  explicit  form  of  (17)  may  be  obtained  straight- 
forwardly but  the  calculations  are  somewhat  tedious  so  only 
the  results  will  be  given  here.  First  lets 

*a  = P Ka  (A-19) 

then  (Op),  (p  = 1,2, 3, 4),  defines  a system  of  coordinates 

on  the  unit  sphere  in  stress  space.  Expressions  for  are  as 
follows : 

^ ^ (1+3  cos  202)  cos  (1-cos  2©2)  cos  2©-^  sin  ©^ 

t, 2 = -^(1-cos  2©2)  cos  2©^  cos  ©^  + i(3 

+cos  2©2)  cos  2©^  cos  2©^  sin  ©^ 

- cos  ©9  sin  2©^  sin  2©^  sin  ©^ 

^ ^ sin  2©2  cos  ©^  cos  ©^  - ^ sin  2©2  cos  ©^  cos  2©^  sin  ©^ 

+ sin  ©2  sin  ©-,  sin  2©^  sin  ©^ 

*4  = ^ U-  cos  2©2)sin  2©-^  cos  ©^  + ^(3 
+cos  2©2)sin  2©^  cos  2©-^  sin  ©^ 

+ cos  ©2  cos  2©^  sin  2©^  sin  ©^ 

Zj  = ^ sin  2 ©2  sin  ©^  cos  ©^  - ^ sin  2©2  sin  ®l  COs  2G3  sin  e4 

- sin  ©2  cos  ©^  sin  2©^  sin  ©^  (A-20) 
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The  ranges  of  the  angles  0 are  as  follows: 

0 < Qx  < 2it  0 < ©^  < n 

0 < 02  < n 0 < 04  < | (A-21) 

The  angles  0 are  suitable  as  a system  of  coordinates  on  the 
initial  yield  surface  of  an  isotropic  material  which  is  given 
by: 


*o  = P<V  <V 


CA-22) 


The  element  of  area  on  this  surface  is  given  by: 
da  = 2Ph  [1  + ( P 1 /P ) 2 1 ^ sin  ©2  sin  3©4  d01  d02  d©3  d©4  (A-23) 


The  unit  normal  vector  to  the  surface  is  given  by: 

[1  + (P’/P)2]^ 


_ P'  Ka 
P 9 64 


(A -24) 


Yield  surfaces. 

In  this  section  a few  well  known  facts  concerning  the 
initial  yield  surface  for  an  isotropic  material  are  restated  in 
terms  of  the  system  of  coordinates  given  in  the  last  section. 

Such  a yield  surface  is  of  course  invariant,  under  those  rotations 
in  stress  space  which  correspond  to  rotations  in  physical  space. 
Obviously  the  yield  boundary  in  the  1,2  plane  determines  the 
whole  surface,  in  fact  the  arc  of  the  boundary  which  lies  in 
the  sector  0 £ ©4  £ ^ is  sufficient  to  determine  the  whole  sur- 
face. Three  of  the  axes  of  the  permissable  rotations  lie  in 
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the  1,2  plane;  one  is  the  x^  axis  and  the  other  two  make  angles 

of  - with  the  x.  axis.  These  three  lines  must  therefore  be  axes 
3 

of  symmetry  of  the  yield  boundary.1 

p _ 

The  yield  surface  must  also  be  convex,  it  is  easy  to 
show  that  if  the  yield  boundary  in  the  1,2  plane  is  convex  then 
this  condition  is  satisfied.  Let  us  assume  that  PC©^)  has  two 
continuous  derivatives.  The  condition  of  convexity  is  explic- 
itly: 


P2  + 2P>2  - PP"  >0  0 < 0^  < | 

P' (0)  < 0 < P» (|) 


(A -2  5) 


The  equality  sign  holds  in  the  first  inequality  for  a straight 
line.  The  Tresca  yield  surface  is  a special  case  of  this  latter 
possibility.  If  the  yield  surface  is  to  be  smooth  the  yield 
boundary  in  the  1,2  plane  must  meet  the  rays  ©^  = 0 and  ©4  = ^ 
at  right  angles,  this  will  be  true  if: 


P = R(cos  36^)  (A-26) 

where  R is  regular  for  0 £ ©^  ^ . If  in  addition  we  are  to 

have  point  symmetry  about  the  origin  then 

P = R(oos2  3©4)  (A-27) 


1 All  members  of  the  sub-group  except  the  identity  have  one 
fixed  axis. 

2 

By  this  we  mean  that  no  part  of  a straight  line  segment  join- 
ing any  two  points  of  the  yield  surface  falls  outside  the 
yield  surface. 
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The  fact  that  the  yield  surfaces  (initial  or  subsequent) 
considered  in  this  paper  are  always  convex  is  intuitively 
obvious  considering  the  fact  that  they  are  all  the  "innermost" 
envelope  of  a set  of  planes. 


Some  additional  spherical  coordinate  systems. 

Several  orthogonal  systems  of  coordinates  on  the  unit 
sphere  in  five  dimensions  may  be  manufactured  as  in  the  follow- 
ing example.  Let  the  projection  of  an  arbitrary  unit  vector  on 
the  (x-^,X2)  subspace  be  of  length  cos  O2  so  that  the  length  of 
the  projection  on  the  orthogonal  complement  (x-^x^jX^,  is 
sin  G2.  Let  the  projections  on  the  x^  and  x2  axes  be: 

= cos  ©2  cos  ei 

(A-28) 


^2  = cos  @2  sin 


as  in  plane  polar  coordinates.  In  the 
projections  on  the  axes  are  made  as  in 
coordinates  to  give: 


(x^jX^jX^)  subspace  the 
ordinary  spherical  polar 


= sin  ©2  sin  ©^  cos  ©^ 

= sin  ©2  cos  ©^  (A-28)' 

Kj  = sin  ©2  sin  ©^  sin  ©^ 


The  ranges  of  the  variables  are: 
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0<61<2ti  0 < < k 

0 <G2  < I 0<Ql+<2ti 


(A-29) 


The  element  of  area  on  the  unit  sphere  is  given  by: 

da  = cos  Q0  sin^  Q~  sin  ©_  dQ,  d©n  d0o  dS,  (A-30) 

^ ^ 3 f 2 j n- 

is  given  by: 

JL  J_  = cos3  Q~  cos  3®t  cos  0_  cos  0O  sin^  0o  (1  - 

P3  3 3/3  2 1 2/3  1 2 2 

- 3 cos  2©^)  + sin2  sin2  0^(sin  0-^  cos  ©2  cos  2©^ 

+ sin  ©2  cos  ©^  sin  2©^)  (a-31) 

Another  system  of  coordinates  is: 


h 

= cos 

61 

^2 

= sin 

el 

sin 

e2 

cos 

e3 

0 

< 

9i 

< 

It 

**3 

= sin 

ei 

cos 

^2 

cos 

©4 

0 

< 

©2 

< 

It 

^4 

= sin 

01 

sin 

«2 

sin 

e3 

0 

< 

G3 

< 

2% 

= sin 

01 

cos 

e2 

s in 

04 

0 

< 

e4 

< 

2k 

da  = 

= sin3 

el 

sin 

©P 

cos 

°2 

d©-^  d©2 

d© 

3 ae4 

3^/2  J = cos  3©n  + 4 sin2  ©i  cos2  9p  [3  cos  3n 

2 p3  3 12  i 1 

+ v/3  sin  ©^  sin  9g  cos  (2©^  - ©^ )]  , (a-32) 
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Appendix  B 


Solving  for  Hff). 


Make  the  xQ  coordinates  dimensionless  by  letting: 


Xa  " £ 3a»ij  °ij 


(B-l) 


where  k is  the  yield  stress  in  pure  shear.  Let  the  yield 
surface  be  p = 1 and  treat  all  planes  equally,  then  the  stress- 
strain  relation  is: 


ea  = 


H(f ) ? da 


(B-2 ) 


The  function  H has  to  be  determined  in  terms  of  experimental 
information,  say  a stress-strain  curve  in  tension.  Let  this 
be  given  by: 


d = Ec  e (B-3) 

X S x J 

where  Eg  is  the  secant  modulus.  The  plastic  part  of  the  strain 
is  ;iven  by: 


e P = 
x 

For  a simple  tension  test  d 
so  (4)  becomes: 


(tr-  ~ w)  d 


E 


E' 


x = /3  X]_  = /I  p 


(B-4) 

and  ©!  = >f3  exp 


ei  = 3(ET  " p = 2F(p)  (b-5) 

In  the  polar  coordinate  system  defined  by  eqs.  (A-32)  the 
distance  of  a plane  from  the  origin  for  the  present  path  is: 
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f = p COS  ©]_ 


l,3 

(B-6) 


(sec  fig.  11)  and  (2)  becomes: 


2F(p) 


P©1  f'J  [>2)CP27t 


o do  Uo 


H(p  cos  ©-^)cos  9-j_sin3©-^sin  ©2Cos  ©2d©4d©2d©2d©^ 

(B-7) 


Ig. 

H(p  cos  0^)  cos  9-^  sin^©-^  d©^ 

Let  r|  = p cos  6-p  then  (8)  becomes: 

pS’(p)  = m2  j H(r) ) r)  (p2  - t]2)  dr] 

This  may  be  solved  for  H by  mere  differentiation  to  give! 


(B-8) 


(3-9) 


H(p)  = [(p^J'/p]'/^2? 


(B-10) 


Reduction  to  loading  lines. 

For  loading  in  the  x-^,X2  plane  it  is  convenient  to 
reformulate  the  problem  in  terms  of  loading  lines  in  that  plane# 
In  the  present  case  of  course,  because  of  symmetry,  exactly  the 
same  analysis  applies  to  say  the  x-^x^  plane  or  any  other  flat 
two  dimensional  subspacc  of  stress  space  which  includes  the 
origin.  The  solution  in  these  other  cases  may  be  obtained  from 
that  in  the  x^,X2  plane  by  a transformation  of  coordinates  or 
even  by  merely  relabeling  the  axes. 

In  the  coordinate  system  (A-28)  (2)  becomes: 


1 
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where: 


°1  = 


e2  " 


it 


f 2u_  _ 

^ H(f)  cos  ©^  d©^ 


o 

2it 


H(f)  sin  ©1  d©1 
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(B-ll) 

(B-12) 


H(f)  =ti  H(f  cos  ©2)  sin2  2©2  d©2  (B-13) 

and  where  f is  the  distance  of  a loading  line  from  the  origin 
in  the  x^,x2  plane.  By  using  (10)  and  afler  considerable  reduc- 
tion (13)  becomes: 

1t_ 

H(p  ) = — I 2 [p  cos  © F'  ( p cos  ©)  + F(p  cos  ©)  ] d©  (B-l4) 

^ Jo 

Formula  (l4)  may  also  be  obtained  directly  using  (11)  and  the 
loading  line  approach  without  resorting  to  (10)  and  then  if  so 
dosired  H may  be  obtained  by  inverting  (13),  in  other  words  it 
is  possible  to  return  to  loading  planes  after  having  determined 
H by  the  use  of  loading  lines.  However  if  one  is  only  interested 
in  loading  in  the  x^,x2  plane  this  is  unnecessary. 

In  terms  of  loading  lines  the  plastic  strain  resulting 
from  a simple  tension  test  is  given  by  (11): 


2F(p)  = 


9 

H(p  cos  ©)  cos  © d© 

0-9 


(B-15) 


where  p cos  9=1,  Figure  11  serves  as  an  illustration  in  this 
case  also  but  v/ith  a different  interpretation  than  before. 
Equation  (15)  may  as  well  be  written: 


_ 

F(p)  H (p  cos  ©)  cos  © d© 

o 


(B-16) 


! f 
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because  H(p)  = 0 for  p < 1,  Equation  (16)  is  not  quite  as 
simple  an  integral  equation  to  solve  as  (8).  (See  [18  ] ) 
However  equations  like  (13)  and  (16)  may  be  solved  by  first 
making  a change  of  variable.  Let 

? p p 

r = p 5 s = p cos  0 


then  (16)  becomes: 


2 r^  F(r^)  = 


^ r H(s^)  ds 


jo  ( r - s ) 2 

This  is  the  familiar  Abel's  eq.  and  the  solution  is: 


(B-17) 


H(rT)  = T s^_F.(s?)  ds 

it  dr  ^ (r  - s)2 


= 2 
it 


r d [s^  F(sh_] 
o (r  - s)~ 


(B-18) 


Transforming  bade  again  gives  (14), 


Broken  line  loading  path. 

The  plastic  strain  resulting  from  any  loading  path  in 
the  x^,Xg  plane  may  now  be  written  down  in  closed  form  in  terms 
of  H,  the  function  H having  been  computed  from  the  stress-strain 
curve  using  (l4).  If  however  one  wished  to  compute  the  strain 
for  a large  number  of  broken  line  loading  paths  it  would  be 
more  economical  to  have  a specific  formula  for  that  case.  This 
formula  has  been  worked  out.  Suppose  the  resulting  yield 
boundary  is  as  shown  in  fig.  12  (whatever  the  loading  path), 
then  the  plastic  strains  are  given  by: 
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el  = F<P0>  + cos(a+P)F  (P].)  “ sign  a [C(a,y)  - F(pQ)] 


- sign  p [C(p,y)  - F(p^)  ] cos  (a+p)  (B-19) 


e2  = sin  (a+p)  F(p^)  - sign  P [C(P,y)  - F(p^)]  sin  (a+p)(B-20) 


where : 


C(a,y)  = -|  sin  |a| 


FCcos  Q/cos  y)  dQ 
1 - cos2acosPo 


(B-21) 


Note: 


lim  G(a,y)  = F(l/cos  y)  (B-22) 

a-^  o 


AT  1 -1  06 


b7 


Appendix  C 


Uniform  expansion. 

In  the  coordinate  system  (A-32)  eq.(31)  with  p = 0 

becomes : 


(C-l) 


(see  Fig. 9).  Carrying  out  the  integrations  on  62>e3’®4  and 
dropping  the  subscript  on  9^  gives: 

r 0 

r9  = ! (p  cos  9 - 1 - r2)  sin^O  d9  ( C —2 ) 

'Jo 

where  p cos  9=1+  Equation  (2)  may  be  written: 

p cos  9 - 1 = 2\k2p  ^ (cos  9 - cos  9)  sin^9  d9  (C-3) 

One  more  integration  gives: 

p cos  9 - 1 = ^ \7t2p  (3-8  cos  9 + 6 cos29  - cos^)  (C-4) 

The  r.h.s.  of  (4)  has  the  same  sign  as  X and  is  a monotone 
function  of  9 for  0 £ 9 £ 2.  as  is  evident  from  (3).  Solving 
(*+)  for  p gives: 

p = [ cos  9 - ^ Xti2 (3  - 8 cos  9 + 6 cos^9  - cos^G?  ] (C-5) 

For  0 < X < co  , p increases  steadily  with  9 and  becomes  infinite 

sometime  before  9 reaches  n/2 . However  for  X < 0 no  matter  how 
small  in  absolute  value  the  denominator  in  (5)  is  always  positive 
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for  0 < G ^ . Thus  p does  not  become  infinite  as  © increases 
from  0 but  reaches  some  maximum  value  and  then  decreases  until 
9 reaches  m/2  at  which  time  the  elastic  region  disappears  al- 
together. The  non-existence  of  solutions  for  arbitrarily  large 
p and  the  possibility  of  two  solutions  for  some  values  of  p is 
untenable,  therefore  negative  values  of  X must  be  excluded.  For 
p fixed  and  X — 7 00  it  is  evident  from  (5)  that  9—^0,  in  other 
words  the  corner  disappears  as  X — 7 co . 

It  is  possible  to  depart  from  a radial  loading  path 
by  turning  the  path  through  a sufficiently  small  angle  without 
unloading  any  planes  which  have  already  been  loaded.  Obviously 
all  planes  which  are  continuing  to  load  must  pass  through  the 
loading  point  and  if  no  [lane  which  has  ever  been  loaded  has 
subsequently  unloaded  then  all  planes  which  have  been  loaded 
pass  through  the  loading  point.  In  this  case  the  solution  is 
that  just  given  for  radial  loading  which  is  unique  for  X J>  0. 

It  follows  that  the  solutions  for  nearly  radial  loading  are  also 
unique.  The  question  of  whether  or  not  there  is  a unique  solu- 
tion for  non-nearly-radial  loading  paths  is  left  open.  The 
permissable  angle  of  turning  for  a nearly  radial  loading  path 
is  determined  in  the  following. 

Suppose  the  loading  has  been  along  the  x^  axis  up  to 
some  point  P and  then  a small  increment  of  load  is  added  to 
reach  the  point  P’  off  the  axis  (see  fig.  13).  As  shown  in 
the  figure,  y has  been  the  angle  of  turning.  The  line  marked  t 
is  the  trace  of  a single  plane  and  this  plane  before  all  others 
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will  begin  to  unload  if  y is  too  large.  In  the  figure  t1  is  the 
trace  of  this  same  plane  after  the  increment  of  load  has  been 
added.  As  shown  this  plane  has  remained  in  contact  with  both 
the  loading  point  and  the  expanding  spherical  portion  of  the 
yield  surface;  therefore  it  has  been  neutrally  loaded  and  y is 
the  maximum  perrnissablo  angle  of  turning.  From  the  figure  half 
the  angle  of  the  corner  at  P1  is  given  in  two  ways: 

S.-Q_de=£-e-da 

2 2 


. * . dQ  = da 

(C-6) 

Also 

dp  = ds  cos  y 

(C-7) 

and 

pda  = ds  sin  y 

(C-8) 

From  (6-8) 

it  follows  cot  y = i 

(C-9) 

P d© 


Using  (5))  after  some  manipulation  (9)  becomes: 

cot  y = p sin  0 + 4 cot  1 © (p  cos  0-1)  ?-f.9.9.s.  Q (C-10) 

^ 3+cos  0 

For  fixed  p as  X— >co,  © — . * . y — 5O  and  the  permissible  angle 
of  turning  vanishes.  As  X-*  0,  p cos  0— » 1 and  cot  y -9  tan  0 
as  it  should.  The  permissable  angle  of  turning  decreases  as  p 
increases  for  fixed  X > 0 and  is  always  less  than  half  the  angle 
of  the  corner.  Of  course  no  part  of  the  nearly  radial  path  need 
be  radial  as  was  assumed  at  the  beginning  of  the  argument,  what 
is  essential  is  that  the  loading  path  proceed  into  the  cone 
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beyond  P with  apex  angle  2y  and  axis  OP. 

Translation. 

A procedure  similar  to  that  in  the  last  section  gives 
p as  a function  of  9 (see  fig.  1,+  ).  The  result  iss 

p cos  9 = 1 + ii,  sin^  ^ 9(3  cos^9  + 9 cos  9 + 8)  (C-ll) 

p 

If  [i  < - 15/i+rt  there  is  the  same  sort  of  objectionable  behavior 
as  in  the  last  section  for  \ < 0.  In  the  present  case  also  it 
is  evident  from  (11)  that  0-+  0 as  [i  -?  oo  for  fixed  p . The 
solutions  for  nearly  radial  loading  paths  will  be  unique  if  it 
can  be  shown  that  there  is  a unique  solution  for  the  case  in 
which  all  planes  that  have  ever  been  loaded  pass  through  a given 
point,  say  P.  When  the  path  is  exactly  radial  and  p > 0 as  in 
fig.  lh-  there  is  not  much  doubt  about  the  uniqueness  of  the 
solution  because  of  the  symmetry  of  the  situation.  But  what  if 
the  loading  path  had  reached  P by  a curve  slightly  above  the  x-^ 
axis  — might  not  the  yield  surface  look  something  like  fig.  15 
while  still  all  planes  that  had  ever  been  loaded  pass  through  P? 
This  can  not  be  case  for  p > 0 according  to  the  following 
argument.  In  the  case  shown  in  fig.  15  the  configuration  of 
loading  planes  is  perfectly  symmetrical  about  the  axis  O'P  and 
thus  so  is  the  distribution  of  direct  motion.  The  displacement 
vector  Aq  discussed  in  the  text  must  be  in  the  direction  O'P  for 
p > 0,  but  this  would  have  displaced  the  spherical  part  of  the 
yield  surface  downward  instead  of  upward  as  it  is  shown  in  the 
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figure.  The  only  possible  position  for  O'  is  on  the  line  OP. 

The  argument  does  not  hold  of  course  if  |x  < 0,  nor  would  it 
hold  for  other  initial  yield  surfaces  than  spherical  ones.  How- 
ever in  the  present  case  it  is  certain  that  the  center  of  the 
spherical  part  of  the  yield  surface  must  lie  on  the  line  OP  if 
all  planes  which  have  over  been  loaded  pass  through  P.  The 
half  apex  angle  of  the  cone  of  total  loading  is  given  by  the 
same  formula  as  (9). 
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Fig.  6.  Limited  path  independence  with  Bouschmger  effect 


Fig.  7 Yield  boundary  after  non  nearly  radial  loading 


